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Abstract 

We derive the kinetic equation for fermions and antifermions interacting with a planar Higgs 
bubble wall during the electroweak phase transition using the 'evenisation' procedure and the 
Wigner formalism for a Lagrangian with the phase of the complex fermion mass rotated away. We 
obtain the energy, velocity and force for the particles in the presence of the Higgs bubble wall. 
Our results using both methods are in agreement. This indicates the robustness of evenisation as 
a method to study quantum corrections to the velocity and force for particles in the Higgs wall 
during the electroweak phase transition. We also derive the transport equations from the zeroth 
and first moment of the kinetic equation. 
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I. INTRODUCTION 



In models of electroweak baryogenesis, baryon asymmetry is created during the elec- 
troweak phase transition at temperatures of around 100 GeV. Since this temperature cor- 
responds to an energy scale that is experimentally accessible, there has been considerable 
interest in ascertaining the details of electroweak baryogenesis models, in particular, the 
kinetic equation or Boltzmann equation for the transport of particles through the Higgs 
bubble wall in a first order electroweak phase transition. 

It is well known that, in the absence of collisions, the characteristics of the Boltzmann 
equation are the single particle trajectory equations, i.e. if /(p, x, t) is a distribution function 
for a system of particles then '^'^^^f^'*'' = 0, where x and p satisfy single particle equations of 
motion. Expanding the total derivative using the chain rule, 

OJ + ^d^f + ^dpf = 0, (1) 
dt dt ^ ^ ^ 

and ^ and ^ can be obtained for the system from the single particle equations of motion 
(Hamilton's equations). However, in relativistic quantum theories it is non-trivial to obtain 
the single particle equations of motion from the Hamiltonian because of interference between 
particle and antiparticle states causing Zitterbewegung. But one may obtain equations of 
motion that have a consistent single particle interpretation by making use of 'evenised' 
operators In Ref. P] we used evenisation to calculate particle trajectories to O(^) in the 
presence of the Higgs bubble wall during the electroweak phase transition with a complex 
mass term given by mjiiljilj+imjip'y^ilj = \m\il)e^^^^ ip. In this article we apply the procedure of 
evenisation to study the interaction of fermions and antifermions with the Higgs bubble wall 
with the phase in the mass term rotated away and with an additional term ^df^Oip'y^'j^ip in 
the Lagrangian. We use this to obtain the energy relation E{p, x) and the kinetic equation. 

Alternatively, in Ref. P], the energy relation, ti;s±(p,x), and a kinetic equation in the 
absence of collisions for particles and antiparticles interacting with the Higgs bubble wall are 
obtained from the equation of motion for the Wigner function G^{x, /c), which is the Wigner 
transform of the Wightman function i (Also see Refs. 0,0].) In this approach, 

particle and antiparticle distribution functions are extracted from the Wigner function using 
a spectral decomposition that separates positive and negative energy states. The particle 
and antiparticle distribution functions are then found to satisfy a coUisionless quantum 
corrected kinetic equation in the presence of the external force field provided by the Higgs 



wall. In Ref. the Lagrangian was as in Ref. P|, i.e. with a complex pseudoscalar mass 
term, as mentioned above. Here we apply the Wigner formalism to the rotated Lagrangian. 
We obtain the energy relation and kinetic equation. These results agree with those obtained 
using evenisation. 

The aim of this article is to establish the robustness of the application of the evenisation 
procedure to study quantum corrections to particle trajectories. Evenisation is a convenient 
technique that allows us to reobtain classical expressions {0{h^) for the velocity and force for 
a relativistic quantum system of particles^]. However, as far as we are aware, the extension 
of evenisation techniques to obtain quantum corrections to the velocity and force for a system 



of particles was first done in Ref. and shown to be consistent with the Wigner formalism. 
Below we work with a different Lagrangian and again find consistent results between the 
two methods. 

The energy relation and kinetics (velocity and force) for a fermion in the bubble wall with 

no 

the real mass term were first derived in Refs. p, 12I| using a WKB ansatz for the components 
of the spinor wave function of the fermion. Refs. p, 0| considered states related with 
chirality while below we consider spin eigenstates. However our results are also compatible 
with the results in Refs. 0, Q] ^• 

One may analyse the interaction of the particles with the background Higgs field with 
either the rotated or the unrotated Lagrangian. The final baryon asymmetry is associated 
with the net axial current generated in the wall and its outward diffusion away from the 
wall. While the Lagrangian is not invariant under the rotation, the axial current {ijj'f^'y^ip) 
is invariant, and so the final asymmetry produced should be the same. 

The energy relation obtained for the unrotated Lagrangian with a complex mass is 

hs 



and can not be easily recast as a mass shell condition. The energy relation for the rotated 
Lagrangian with a real mass is 

hsO' 



^ = Vp/ + m2-^ (3) 
and can be more easily restated as a mass shell condition (E + _ ^2 _ This is 



^ Eq. (9) of Ref. |3 is missing an overall minus sign in the expression for the force in the z direction. 
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similar to the mass shell condition for classical systems and may make it easier to interpret 
the behaviour of the particles of the system. 



II. KINETIC EQUATION USING EVENISATION 



A first order electroweak phase transition proceeds via the formation of Higgs bubbles. 
As the bubbles expand they move through the ambient sea of quarks, leptons and other 
particles. In this section we first obtain a semi-classical expression for the single particle 
energy in terms of the evenised position and momentum. We then obtain expressions for the 
evenised velocity and force and substitute these in Eq. (jT)) to obtain the kinetic equation. A 
brief introduction to evenisation is given in Ref. Here we merely present the necessary 
formulae for our analysis. 

Any operator A can be split into an even part A and into an odd part {^}- Even and 
odd parts of the operators are defined by using the sign operator 



A 



H 



(4) 



where H is the Hamiltonian. The eigenvalues of the sign operator are ±1, corresponding to 
particle and antiparticle states. Then 



A 



^ (i + aAa 



{A} = 1{A- aAa 



The even part of the product of two operators A and B can be written as 



AB 



A 



B 



(5) 
(6) 

(7) 



We now apply the above to the system under consideration. 

The Lagrangian describing the interaction of particles with the bubble wall can be mod- 
eled by 

— 1 — 771 — 

c = iii)^^ + -d^ei^j^j^ij - (8) 

This is a rotated form of the Lagrangian of Ref. j2| with m above equivalent to |m| in Ref. 
[2]. The Higgs bubble will be treated as a background field which provides in the bubble 
wall frame a spatially varying real mass for the particles, and a term associated with the 



axial current. We shall consider the limit of large bubbles when the walls can be treated as 
planar. The corresponding Dirac equation 

1 7T7 

{^0+-^^erl'-J-)i^ = Q (9) 



can be rewritten in the form 



ih^-^ = , (10) 



with the Hamiltonian H given by 

H = a^if + ^d^ef) + /3m (11) 
= dy + ^mS' + /3m . (12) 

We have defined the spin operator S as 07^, without the standard h/2 in the definition 
to be consistent with the notation in Ref. j^; the spin operator defined here is S with 
Pauli matrices along the diagonal. We have adopted the metric (1,-1,-1,-1). In this section' 
designates an operator. We now use the method of evenisation and evenise H^. We discuss 
later why we evenise rather than H. 
We first write after some algebra, 

= p2 +^2 (QiQsiY^g^ip p^^^ +nd^emS'j3 + -a'S^\p\dW] + h&eS^a'p' (13) 

We work in the wall frame and take the wall to be planar in the x — y plane and so m is a 
function of z only. Then, keeping only terms to 0{h), the Hamiltonian can be written as 

1/2 = p2 + ^2 _ s^^f^^ind^m) - he'ms^p - ne'p's^a' , (i4) 

where z = x^, dz = d/dx^ and d' = dzO = —d'^d. We chose to express a^/3 in terms of as 
this simplifies the evenisation of the corresponding term below. 

We now evenise [^2] to order h. The spin operator in general does not commute with 

n 

the Hamiltonian. In order to simplify the problem, as in Ref.|3|, we assume that we are 
in an inertial reference frame where x and y components of the momentum are zero and 
consequently and will be absent from the Hamiltonian. This will allow us to set {S^} 
to later. As mentioned earlier, we wish to express [H'^] in terms of [pz] and [z\ (where 
Pz = p^). The odd part of Pz will be proportional to [A,!)^] (see Eq. © above) and so will 
be of order h. Therefore, using Eq. we can approximate [pi] by [pz]^- Similarly we 



approximate [m(z)^] as [m(i)]^. Furthermore, since {z} is 0{h), ['m{z)] = Tn{[z]) + 0(^^) 
as any dependence of the even operator [^(-2)] on the odd operator {z} can only appear as 
{z}"^- (Expand m(z) in a Taylor series about z = and then replace z hj z = [z] + {z}.) 
The last three terms on the right hand side of Eq. p4|) are 0{h) and therfore we will use 
the sign operator defined upto 0{fi^) to evenise these terms. We define the zeroth order 
Hamiltonian and energy as 

= a^p^ + (3m (15) 



^0 = vpTT^ (16) 

where Pz'^ and are real numbers and are expectation values of the corresponding operators 
in an eigenstate of definite energy and spin. With this we define the sign operator upto 
O (h^^ as follows 

^0 = -7=4=^ (17) 

Note that (Aq)^ = l + 0{?i) which we shall use in the derivation of evenised operators in the 
Appendix. We obtain 

[ff] = 0{h) (18) 
[/3] = ^Ao + Oih) . (19) 

= I^Ao + 0{h) . (20) 

Now 

[S^] = [a^f] = , (21) 

i.e., {5*^} = as S*^ commutes with the Hamiltonian (which is true only in the chosen inertial 
frame). Then, keeping in mind that the odd parts of m, dztfi, 9' and p^ that appear in the 
last three terms of are 0{h), we get, to 0{h), 

[H'] = [pz]' + m' - M'^S'ko - M'M^^^Ao . (22) 

where m and 9' are now functions of [z\. Replacing evenised operators [z\ and [pz] by real 
numbers representing their expectation values in a state of definite energy and spin, \E,s) 
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for particles and | — E, —s) for antiparticles ^ , we deduce the corresponding expression for 
the energy to be 

^2 ^ p^2 ^ ^2 _ f^gQ'fJ^^ (23) 

where we have also replaced by its eigenvalue ±s and Aq by ±1 for particle/antiparticle 
states. The energy relation for particles and antiparticles respectively can then be written 
upto order h as 

hsO' 



E = ^p,2 + ^2 _ __ (24) 

Note that if we had wished to obtain the energy relation by evenising H instead of evenizing 
we would have faced a problem as we would have needed A to 0{h) which itself requires 
E to 0{h). However, since {H} = 0, [H] = [H^Y^^. We are now able to define the sign 
operator upto 0{h) as A = H/E and can use the same to obtain the evenised velocity and 
force to 0{h) for the kinetic equation. 

We now obtain the kinetic equation to 0{h). Using the chain rule for partial differentia- 
tion the kinetic equation is written as 



dtfs± + ^djs± + ^d,Js± = , (25) 



where fs± are the particle and antiparticle distribution functions. We will associate dz/dt 
and dpz/dt with the expectation values of [dz/dt] and [dpz/dt] in states of definite energy and 
spin, as before. Implicitly we are assuming here that the form of the quantum Boltzmann 
equation is the same as that of the classical Boltzmann equation and that the quantum 
corrections are contained only in the coefficients of the equation, namely, in the expressions 
for the velocity and the force. The evenised expressions for the velocity and force to 0{h) 
are derived in the Appendix as 



[dz/dt] = [-^[z,H]' 



E 2El 



A (26) 



Jdt] = [-^[p.,H] 



^ In the earlier version of this paper, we had identified antiparticles of spin s with negative energy particle 



tipa 

solutions of spin s, rather than —s (see Sec. 7.1 of Ref. [ij). Hence the expressions for energy and force 
were different for particles and antiparticles of the same spin. Our current results reflect the P and CP 
violation and C conservation properties of the Lagrangian in Eq. ||HJ). 
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+ — ^ A. (27) 



\ 2E iEi 2 ^ 

Substituting the eigenvalues of the above operators in the kinetic equation ^, we get 

m^' / ne's\ hsO"' 



1 



dpJs± = Q. (28) 



2E \ 2EoJ 2 

Simphfying the above equation by using Eq. to express E as Eo[l — TiO's/ (2£'o)] we get 

Several comments are in order here. If we apply Hamilton's equations to the energy 
relation in Eq. (j21I), i-e., dz/dt = dE/dpz and dpz/dt = —dE/dz, we get expressions 
for dz/dt and dpz/dt that agree with the expectation values of the corresponding evenised 
operators [dz/dt] and [dpz/dt]. But the real number Pz appearing above in Eq. ()24j) need not 
be the canonical momentum. It represents the expectation value of [pz] whereas the canonical 
momentum would be associated with the expectation value of Pz- Secondly, ignoring for 
now the distinction between pz and [pz], Ehrenfest's theorem for a relativistic system would 
imply d{z)/dt = {dH/dpz) and d{pz)/dt = —{dH/dz), and not d{z)/dt = d{H)/d{pz) and 
d{pz)/dt = —d{H)/d{z). Nevertheless the energy relation in Eq. ()24|) behaves as a quantum 
corrected classical Hamiltonian, i.e., one may apply Hamilton's equations to it to obtain the 
velocity and the force. However, in general, one can not presume that E is the quantum 
corrected classical Hamiltonian and indeed in Ref. P| with the unrotated Lagrangian one can 
not treat E{pz,z) as the quantum corrected classical Hamiltonian with pz as the canonical 
momentum. 

In the literature 0, 0] there has been discussion on identifying the canonical and kinetic 
momentum of the particle. One may presume that since Pz is identified with the —idz 
operator it is conjugate to z and represents the canonical momentum. Furthermore pz 
acts as a canonical momentum if the energy relation of Eq. (f^^ is treated as a classical 
Hamiltonian. Now our energy relation also agrees with with that obtained below using the 
Wigner formalism. But that energy relation is necessarily in terms of the kinetic momentum 
as the Wigner function is always written in terms of the kinetic momentum of the system (see 



^ The overall minus sign in the expectation value for the velocity and force for antiparticles can be absorbed 
by taking the expectation value of the momentum in antiparticle states to be —pz, as for systems with 
constant mass, where the energy eignenstates are also momentum eigenstates. 
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Sec. 3 of Ref. P]). This indicates that our real number pz = {[pz]) in the energy relation 
above is the kinetic momentum of the system. Thus the evenised canonical momentum 
in this study becomes equivalent to the canonical momentum and the kinetic momentum, 
unlike in Ref. 0]. 

The single particle trajectory has been obtained in a frame in which the x and y compo- 
nents of the momentum are 0. For a single particle such a frame may be defined but not 
so for a collection of particles. Therefore the treatment above may be treated as that for a 
1+1 dimensional universe only. It is non-trivial to derive as above the kinetic equation for 
a 3+1 dimensional universe because of the presence of spin. 

III. KINETIC EQUATION USING THE WIGNER FORMALISM 

We now present a more formal field theoretic derivation of the kinetic equation in a 
manner analogous to that adopted in Ref. for the unrotated Lagrangian. As stated in 
the Introduction, one can obtain the kinetic equation from the equation of motion for the 
Wigner transform of the Wightman function. We shall closely follow the derivation in Ref. 
0] for the unrotated Lagrangian. The Wightman function is defined as 

G<p{u,v)^z{Mv)ijM) . (30) 

Starting from the Dirac equation for the fermions in the electroweak plasma, one can obtain 
the equation of motion for the Wightman function by multiplying the Dirac equation by iip, 
passing it through the differential operator and taking the expectation value. 

[^ K + ^7V]p.^ (,Uv)^M) = (31) 

Subsequently, performing a Wigner transform of G^{u, v), one obtains G^{x, k). 

g<{x,k) = J d^re^''^>'G<{x + r/2,x-r/2), (32) 

where r = u — v. Unlike the Wightman function which is a function of two spacetime 
variables, Q'^{x,k) is in a mixed representation, i.e., it is a function of x and k and hence 
can finally give a kinetic equation for particle and antiparticle distribution functions which 
are functions of position and momentum. Q^{x, k) satisfies 

rh{x) + ^7^7'1 = 0, (33) 







ih 




= + 




171 






die 


= id,e] 





where 

(34) 

(35) 
(36) 

(Clearly, " in this section is defined differently than in the previous section.) We now assume 
that the bubbles are large enough that we can take the wall to be planar (and in the x — y 
plane). Then, working in a wall frame, as in Ref. P|, in which the momentum of the particles 
parallel to the wall vanishes, S^i'j^Q'^ = si'-f^Q^. Therefore, since (/ — sS^)i'~f^Q^ = 0, 
G^{x, k) may be expressed, in the chiral representation, as 

-^l'Q^ = gf^lil + sa') (37) 



^[9oI + 9tp')^^iI + scr'), (38) 

where a* = p' are the Pauli matrices and Qq^ are functions of k^,x^,t and s. ^ As discussed 
earlier, the Wigner function above is now for a 1+1 dimensional universe. 

This decomposition allows one to reduce the 4x4 equation for the Wigner function, Eq. 
to a 2x2 equation for of. Expanding the operator acting on —i'-f^Q^ as a direct product 



of 2x2 matrices (see Ref. 



we get 



[p + ^qIq _ ^pp3 _ ^pi^g< ^ . (39) 

(Here we have used the definition of the gamma matrices as in Refs. P, 0], which differs 
from the definition of 7" and 7^ in Ref. [lO| bv a minus sign. The metric is (1,-1,-1,-1).) 



The corresponding equation obtained in Ref. [2| for the unrotated Lagrangian is 

[P - sPp^ - mop^ - rh5p^]gf = (40) 

where 

mo,5 = mR,7e-*^^-^\ (41) 



Here we have used the definition of the direct product as in Ref. jj, [lOj . The order of the two terms in 
the product is reversed in Ref. 1^. Also, the above decomposition of G'^{x, k) is as in Refs. which 
differs by a factor of 2 from the decomposition in Ref. 



9 



where mn and m/ are |m| cos^ and |m| sin^ respectively. Expanding gf in terms of I and 
the Pauh matrices as in Eq. ipH)) and muhiplying Eq. ()39|) by I and p* and taking the trace 
gives 







sPg'^- mg{ 


= 


(42) 






sk^gl - img^2 


= 


(43) 


(P + 


+ 


isk^g-2 - mg-Q 


= 


(44) 


(P + 




isk^gl + irhg^ 


= 0. 


(45) 



These equations are similar to Eqs. (2.16)-(2.19) of Ref. ^ but with A;° ^ A;° + fdle, 
rriQ — > m and rh^ —>■ 0. One can independently set the real and imaginary parts of these 
equations to 0. The imaginary parts contain time derivatives and so provide the kinetic 
equations. The real parts provide constraint equations. Assuming that m and d^O vary 
slowly in the Higgs bubble wall we may expand m and 536* to 0{h) as 

ih 

rh = m + —m'dk^ (46) 

die = e' + 'h"d,^ . (47) 

Henceforth, we define z = and kz = k^, and / denotes d/dx^ = d/dz and d^^ = d/dk^. 
dx ■ dk = — dz ■ dk^. The constraint equations then are 

(fc° + y ^')^7o - skzgl - mgl = (48) 
(fcO ^ _ ^f^^gs ^ ^^'d^y^ = (49) 

{k' + 'p')gl + '^dzg', -mg^, = (50) 
ik' + y^')^72 - Y^^9l - \ni'd,^gl = 0. (51) 

The above equations imply relations between gQ and g^. Using Eqs. (j49l5(JI51|) iteratively 
to 0{h) we express the g^ in terms of g^. Substituting this in Eq. (jlHj) gives 

[{k' + ^er-k'z-m^^ 9-0 = 0,- (52) 

This implies 

Ql ^ {k' + ^e'f - kl -m'=0, (53) 
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whose solutions are = ±uJs±{z^ kz) with 

fis I 

uJs±{z, k^) =UqT —d' , ujQ = \Jkl + m? (54) 

uJs+{z, kz) and c<j(_s)_(2;, —k-^) will represent the energy of particles/antiparticles as they pass 
through the Higgs bubble wall. This agrees with Eq. (j24j) . 
The kinetic equation for is 

dtg', + ^e"d,.y, + sdzQl - m'd,^g{ = . (55) 

Using the constraint equations to re-express in terms of g^, the kinetic equation for g^ 
implies 

^^^0 + kA [-^^) - Im^'d.. [-^^) + Y^"^^.^o = . (56) 
Analogous to Ref. Eq. implies that gQ may be written as 

g^, = AnNlk^ml) . (57) 

However, for reasons discussed below, we define 

g'o = AnNlk^ + hse'/2\S{Q^) (58) 

.^4.7V^^±5£^.(.«T..J. (59) 

± z,^o 

is a function of k^, kz, z, t and spin s. 

One now substitutes for gQ from Eq. in the kinetic equation and integrates over 

positive and negative energies separately to get 

s,/..+ g)a./,.+ (-g + ^^r)a,/.. = o. (eo) 

where fs±ikz,z,t) are functions given by 

fs+ = N{uJs+,kz,z,t) (61) 
fs- = 1- N{-iu^_s)-,-kz,z,t) (62) 

If we can identify /s± with the particle and antiparticle distribution functions ^, the above 
equations have the form of the classical kinetic equation, with quantum corrections in the 
coefficients of the derivatives of fs±. 



See footnote Pi- 



ll 



The identification of tlie particle and antiparticle distribution functions can be made in 
two ways. One can argue that the distribution function is that quantity which gives the 
number density on integration over momenta. Defining the current as 

fix) = {ij{x)rij{x)) (63) 

we rewrite it as in Sec. 4.1 of Ref. as 

fix) = -Tr[7^G'<(x,a:)] = - / ^^IY^G^ik,x) (64) 



Substituting for the Wigner function from Eq. (jSH)), and using TiA ^ B = TtA TtB and the 
tracelessness of the Pauh matrices, we get 

f = ^J-^.9oik',k.,z,t). (65) 
If we use the decomposition in Eq. for we find that we can write j° as j° — where 

j1= T.fs±= E (66) 

s=±l s=±l"' l^^J 

(We have ignored the vacuum contribution to j^.) Thus if we define particle and antiparticle 
number densities (per unit length) to be it seems appropriate to identify fs± with particle 
and antiparticle distribution functions. Alternatively, one can argue that the quantum 
corrected kinetic equation should have the same form as the classical kinetic equation, i.e., 
it should contain only derivatives of the distribution function and not terms proportional to 
the distribution function. Since the equation for fs± satisfies this criterion it again seems 
correct to identify fs± with particle and antiparticle distribution functions. However, note 
that if we had used Eq. (jKTj) as the decomposition for neither of the above criteria would 
be satisfied. This is why we chose the decomposition in Eq. (j59p . We further note that with 
the fs± identified as particle and antiparticle distribution functions Eq. (jUnjl above agrees 
with the kinetic equation obtained in Sec. II. 

We briefly comment on the application of the above ideas to the exposition in Refs. 
j^, y, 0] for the unrotated Lagrangian. While fs+ in these works satisfies the second criterion 
above, i.e., the quantum corrected kinetic equation has the form of the classical kinetic 
equation, the quantity whose integral over momentum gives j° are i.e.. 



Js+ 



/^/.+/^.+ , (67) 
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where 




1 



(68) 



s+ = 



(See Eq. (2.41)ofRef. |3|] but with the normalisation of (7q as in Ref. ^) If one rewrites 
the kinetic equation in terms of the variable fs+/Zs+ one gets an additional term in the 
kinetic equation. 

But let us now consider the quantity which gives the number density of particles on 
integration over 4- momentum, or rather, 2-momentum in our 1+1 dimensional case, after 
multiplying with kinetic energy and a delta function containing the mass shell condition, 
i.e.. 



where is the kinetic energy. When converting 5{VL^) to 5{k^ — ujs+) one ultimately gets 
a factor of Zg^ in the denominator. (For the unrotated Lagrangian, the kinetic energy 
equals the total energy ujgj^, as evidenced by the definition of the kinetic momentum in Eq. 
(2.14) of Ref. as uJVg, where uj is the total energy and Vg is the velocity.) Now if we choose 
to define the particle distribtuion functions as -Fs(x, k) with replaced by ujg^, i.e., in terms 
of the function in the integrand of when expressed as an integral over k^ and as in 
Eq. rather than as an integral over only kz as in Eq. (jEZj), then we can identify fs+, 

rather than with the particle distribution function and one also has a 'standard' 

form for the kinetic equation for this particle distribution function. Therefore this prompts 
us to modify our first scheme of identification of the particle distribution function. It does 
not affect the case of the rotated Lagrangian. 

The transport equations relevant for electroweak baryogenesis for the unrotated La- 
grangian are obtained in Eqs. (5.2), and (5.8) and (5.9) of Ref. \^ from the kinetic equations 
for Qq and g^, which agree with taking the zeroth and first moments of the kinetic equation 
for fs^, i.e., by integrating the kinetic equation over kz after multiplying by 1 and the veloc- 
ity V. For the rotated Lagrangian considered in this article we obtain the kinetic equation 
using both methods and they agree (for us, v = kz/uj^). The equations we obtain are 



^ The Z factor is absent in Eq. (4.10) of Ref. |5j because the quantity 5fg in the numerator is aheady 




(69) 



dtns± 



+ dz{ns±Us±) = 



(70) 



o{n). 
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dt (n,±M,±) + (n,±(t;2^)) + -m^^.^ - -s{m^e')%s± = (71) 

where 

f dk, 
= J 27^^^ 

and Us± = {vs±), and 

It is interesting to note that the form of the above transport equations for the rotated 
Lagrangian are identical to those in Ref. jy for the unrotated Lagrangian (if one amends 
Ref. ly to identify antiparticles of spin s with negative energy particles of spin —s). 

IV. CONCLUSION 

In conclusion, we have obtained the energy relation and the kinetic equation for fermions 
and antifermions interacting with the Higgs bubble wall during the electroweak phase transi- 
tion using the method of evenisation and the Wigner formalism. Our results for the velocity 
and force on the particle/antiparticle as they pass through the Higgs wall are the same using 
both methods. This indicates that evenisation is indeed a reliable method to investigate the 
quantum corrections to the velocity and force acting on particles and antiparticles as they 
traverse the bubble wall during the electroweak phase transition. 
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Appendix 

We present below the derivations of the various evenised expressions quoted in the text. 
Standard formulae that we shall use are (3"^ = — (7°7^)^ = (a^)^ = 1. The symbols square 
brackets are used below for both commutators and evenised operators but the difference is 
clear from the context. 

Evenised (3, and 7^7^ to 0{h°): 
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(72) 

fs±, (73) 
= 2,3). (74) 



To o(h'^), 



(3 + Ao/9A, 







Eq 

= T^Ao, (75) 

where we have used Aq = 1 + O(^) in the last equahty. Since the Ihs is even and the sign 
operator on the rhs is even, m.{z) may be replaced by [m(i)], which, as explained in the 
text, is equal to to 0{h). Therefore 

[/3] = ^Ao + 0(^). (76) 



Similarly, to 0{h^), 



2[a^] = + Aod^Ao 



a-^ H ^——^—a'^Ao 



Eq 

,3 

Eq 



a-^ + a-^ ' — Ao 



^ f A„ ^ (77) 

Again, since the Ihs is even and the sign operator on the rhs is even, may be replaced by 
\pz\- Therefore 

:,3i _ \Pz^ ■ 



Analogously, 



Therefore, 



[a^] = ^Ao + 0(^). (78) 



2[7°f] = ff + KffK 

= 0. (79) 

[ff] = 0{h) . (80) 
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Evenised expressions for the velocity and the force to 0{h): 
Since we wish to work till 0{h) we now use the sign operator A defined as 

A 4 (81) 

with E defined to 0(h) as in Eq. (j21|). Now dz/dt = -{i/h)[z,H] = dr". Then 

2[a^] = a^ + Aa^A 

= a^ + a-^ ^- A 

= ^ ^ = (^^^ 

Since the Ihs of the above equation and the sign operator on the rhs are even we may 
replace by [pz]- Since is also even, we may similarly replace by [^'d^]. Since {9'} 
is proportional to [^', A] ~ we may rewrite as [^'][d;^]. As discussed in the text for 

[§'] = 9'{[z]) + 0{h^). Further we may use the expression for [d^] obtained above in 
Eq. (f7S|) upto 0{h^) in the rhs of the above equation. Therefore we get, to 0{h), 

ld^/dt]=[-^- JA. (83) 

The force dpjdt = -{i/h)[pz, H] = -f3m' + h9"S^/2. Then 

2 [/3m'] = {3m + Ai3m A 

= pm + pm A — pa — A 

^ ^ E ^ E 

2m'm — hf3m'9'S^ + ih^a^rh'' ^ , , 

= A. (84) 

Again, we may replace m'm by [m'm] as both sides of the equation should be even. As 
{m'}{m} is 0{h'^), Eq. ((Tj) implies that this reduces to [m'(£)][m(5)]. This may further 
be rewritten as m'([z])m([i]). For the second term in the numerator, we replace j3rh'9'S^ 
by [l3rh'9'S'^\ = [Pm'9']S^. This may be rewritten as [/^m'] [^'jS*^. As above, we replace 
[^'(i)] by 0'([i]). Since, by Eq. (jH)), {rh'} is 0{h), [Pm'] can be rewritten as [/3][m'], and 
[P] may be replaced by the 0{h^) expression in Eq. (f7H|) and [m'] by m'([z]). For the third 
term, we replace (3a^m" by [(3a^m"]. Since this term is already of 0{h) we can write 
it as [/3Q;^][m(i)"] and further as [/3d^]m"([i]). Now (3a^ = S'^7°7^ and since {S'^} = 0, 
= 5'^[7°7^] = 0{h). Therefore, this term is of 0{fi^) and we may ignore it. Thus 

l^,-„T. W_M1!!»M!)a. (85) 
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Therefore 



m 



2' nm^'e'A.s'^\ , hff's^ 



2E ' 4:El 



m and 6 above are functions of [z] 
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